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Detection  of  single  or  multiple  optical  subsurface  laminar  features  (e.g.,  thin  layers)  in  marine  waters 
has  many  implications  on  ecological  studies,  management  of  fisheries,  and  military  applications.  This 
study  has  four  objectives:  1)  to  corroborate  a  previous  index  based  on  remote  sensing  reflectance  ratio 
(Rl=  Rrs(443)/Rrs(490))  for  predicting  depth  of  lidar-derived  baekseattering  layers,  2)  to  evaluate  at 
what  extent  these  relationships  hold  in  different  marine  regions,  3)  to  examine  the  nature  of  R1 
variability  in  terms  of  inherent  optical  properties,  and  4)  to  investigate  the  influence  of  water 
stratification  on  spatial  distribution  of  R1  values.  Measurements  of  inherent  optical  properties 
(absorption  eoeffieient  and  scattering  coefficient)  were  obtained  from  a  towed  underwater  vehicle 
(Scanfish)  in  three  geographic  locations,  two  in  US  (Monterrey  Bay,  MB,  and  East  Sound,  ES)  and  one 
in  Turkish  (Black  Sea,  BS)  coastal  waters.  For  each  site,  ease  studies  were  examined  based  on 
subsurface  optical  layers  distributed  at  two  different  depths  (ES:  7  and  20  m,  MB:  7  and  20  m,  and  ES: 
15  and  27  m).  R1  was  theoretically  derived  from  each  Seanfish  profile  and  the  R1  skewness  (v|/)  was 
calculated  for  each  case  study.  The  magnitude  of  the  total  absorption  eoeffieient  at  675  nm  suggested 
large  differences  (>100%)  in  trophic  status  between  the  studied  areas.  However,  a  common  statistical 
pattern  consisting  of  lower  \|/  -  deeper  optical  layer  was  found  in  all  study  eases.  This  variation  was 
explained  by  optical  differences  above  the  ‘optieline’  and  mainly  related  to  changes  on  scattering 
eoeffieient  of  particulates.  In  general,  skewness  of  R1  was  more  influenced  by  ^>(440)/o(440)  rather 
than  by  a(488)/6(488)  spatial  variability.  The  observed  relationships  between  v|/  and  the  depth  of  the 
optical  marine  layer  confirmed  a  previous  finding  using  airborne  lidar  and  ocean  color  data.  Also,  our 
analysis  supports  the  use  of  V|/R1  to  identify  the  relative  vertical  position  of  optical  submerged 
structures  in  marine  environments  with  variable  coupling  between  water  stratification  and 
phytoplankton  layers. 

INTRODUCTION 

A  common  assumption  of  remote  sensing  algorithms  based  on  ocean  color  sensors  is  the  vertical 
homogeneity  of  the  water  column  in  terms  of  optical  properties.  This  approximation  is  very  often  not 
met  in  coastal  and  oceanic  stratified  waters  due  to  the  presence  of  laminar  features  altering  the 
underwater  light  field.  These  submerged  layers  commonly  correspond  with  thin  layers  (i.e.,  <3  m 
thick),  have  typically  high  concentration  of  dissolved  and  particulate  material  with  respect  to  the 
surrounding  medium,  and  arc  preferentially  developed  along  the  horizontal  component  (Chumside  and 
Donaghay,  2009).  Technologies  used  to  detect  vertical  location  of  subsurface  optical  layers  are 
commonly  based  on  lidar  (Light  and  Range  detection)  systems  (Hoge  et  ai,  1988;  Barbini  et  al.,  2003; 
Chumside  and  Donaghay,  2009).  Unlike  these  investigations,  in  a  recent  contribution,  we  showed  a 
new  approach  to  discriminate  waters  with  shallow  versus  deep  optical  layers  based  on  passive  optical 
measurements  (Montes-Hugo  et  al.,  2010).  Briefly,  the  relative  distance  of  the  optical  layer  to  the  sea 
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surface  is  estimated  by  calculating  the  third  moment  around  the  mean  (i.e.,  skewness  or  \j/)  of  a  specific 
remote  sensing  refleetanee  ratio  (Rl=  Rre(443)/Rrs(490)).  As  v|/Rl  decreases  the  subsurface  optical 
layer  becomes  deeper.  The  present  study  has  four  main  objectives.  First,  it  will  test  optical 
relationships  found  by  Montes-Hugo  et  al.  (2010)  using  an  alternative  approach  consisting  in 
independent  calculations  of  \|/R1  based  on  in  situ  measurements  of  inherent  optical  properties  (lOPs) 
obtained  with  an  undulating  ROTV  (remotely  operated  towed  vehicle).  Second,  it  will  examine 
whether  the  aforementioned  approach  ean  be  generalized  across  different  marine  coastal  domains 
located  at  different  latitudes  and  eharaeterized  by  different  trophic  status.  Third,  it  will  investigate  the 
mechanisms  explaining  \|/R1  changes  in  terms  of  lOPs  modifications.  Lastly,  it  will  quantify  the 
influence  of  water  stratification  on  biological  ‘layering’  and  its  impact  on  \|/R1  patterns. 

METHODS 

ROTV  surveys 

High  vertical  resolution  profiles  (0.02  to  0.56  m)  of  in  situ  optical  (i.e.,  total  absorption  eoefficient  of 
dissolved  +  particulate  matter,  a,  and  beam  attenuation  eoeffieients,  c)  and  CTD  measurements  were 
obtained  using  a  ROTV  (Seanfish  MK II,  intelligent  undulating  vehicle,  EIVA,  Denmark)  in  coastal 
oligotrophie  (Black  Sea,  BS,  41.22-43.70  °N,  28.90-31.26°E),  mesotrophie  (Monterrey  Bay,  MB, 
36.33-36.83  °N,  121.05-122.85°W)  and  eutrophic  (East  Sound,  ES,  48.62-48.67  °N,  122.86- 
122.89°W)  waters.  The  ROTV  (0.9  x  0.3  x  1.8  m,  length,  height,  wide)  was  operated  using  an 
undulating  mode  complemented  with  a  wineh  to  expand  horizontal  sampling  range  up  to  400  m.  The 
ROTV  has  a  weight  of  75  kg,  a  maximum  depth  range  of  400  m,  and  a  payload  of  approximately  50 
Kg.  ROTV  data  were  collected  between  14  and  16  pm  (local  time)  and  two  case  studies  were  selected 
per  study  site:  ‘shallow’  (NS)  and  ‘deep’  (FS)  subsurface  optical  layer.  The  ROTV  diving  ‘saw  tooth’ 
pattern  varied  between  studied  areas  (Fig.  1).  Depth  interval,  vertical  resolution  and  horizontal  spacing 
between  profiles  were  2.8-65,  0.25-0.56  and  370.9  m,  for  BS,  1.6-48,  0.36  and  268.5  m  for  MB,  and 
0.7-22,0.02  and  87.4  m  for  ES,  respectively. 


Horizontal  distance  (m) 

Figure  1.  An  example  of  Seanfish  MK  II  dives  (ups  and  downs)  in  the  Black  Sea  j^BS),  Monterrey 
Bay  (MB),  and  East  Sound  (ES).  BS,  MB  and  ES  sampling  locations  had  an  averaged  bottom  depth 

of 125, 120  and  27  m,  respectively. 
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Analysis  of  optical  and  CTD  profiles 

Raw  a's  and  c’s  coefficients  were  derived  from  an  ae-9  (WetLabs)  at  three  wavelengths  (440,  488  and 
675  nm)  and  corrected  by  salinity  and  temperature  effects  (Pegau  et  ai,  1997;  Langford  et  ai,  2001) 
using  CTD  (SeaBird  911)  measurements.  Scattering  residuals  were  removed  (Zanaveld  et  ai,  1994). 
Only  descending  dives  including  ROTV-derived  lOPs  (i.e.,  total  absorption  and  beam  attenuation 
coefficients),  temperature  and  conductivity  were  processed  and  smoothed  every  1-m  along  the  vertical. 
Salinity  and  seawater  density  at  each  depth  was  derived  from  CTD  variables  (i.e.,  temperature  and 
conductivity  without  pressure  correction)  and  using  the  standard  UNESCO  polynomial  equation  of 
state  (Millero  et  ai,  1980;  Fofonoff  and  Millard,  1983).  For  each  variable,  data  gaps  were  removed  by 
linear  interpolation.  Missing  determinations  were  more  common  at  finer  spatial  sampling  frequencies 
(e.g.,  in  ES).  The  upper  mixed  layer  depth  (UML)  was  computed  based  on  a  density  difference 
threshold  with  depth  of  0.015  Kg  m"^. 

Modeling  of  ocean  color 

Spectral  remote  sensing  reflectance  (Rrs(A.))  was  estimated  for  each  descending  profile  based  on 
vertical  distribution  of  a  and  c  values  at  specific  wavelengths.  A  light  propagation  model  (Hydrolight, 
Sequoia  Inc.)  and  ancillary  data  provided  from  airport-based  meteorological  stations  (Black  Sea, 
www.infospacc.ru,  Monterrey  and  Eastsound,  www.wunderground.eom)  were  used  to  simulate 
Rre(440)  and  Rre(488).  Spatial  skewness  (\|/)  in  each  experiment  was  calculated  for  theoretical 
Rre(440)/Rre(488)  ratios,  and  vertically-averaged  optical  properties  measured  above  the  subsurface 
layer  (e.g.,  a(675)).  Contribution  of  total  scattering  (b  =  c  -  a)  with  respect  to  absorption  coefficient  to 
V|/R1  was  also  investigated.  First  optical  depth  (Zoo)  was  calculated  as  the  inverse  of  Kd(488)  or  the 
vertical  attenuation  coefficient  of  downwclling  diffuse  light,  with  Kd(488)  ~  a(488)/jy  (Mobley,  1994), 
where  ^  is  the  average  cosine  and  was  equal  to  1  based  on  cloud  cover  conditions. 

RESULTS  AND  DISCUSSION 

Cross-sections  of  a(675)  suggested  that  phytoplankton  was  an  important  optical  component  in  all 
subsurface  layers  under  investigation  (Fig.  2).  Consistent  with  an  increase  of  solar  radiation  and 
freshwater  river  discharge  as  the  spring-summer  season  progresses,  the  main  a(675)  layer  and 
pycnocline  were  always  deeper  during  the  second  experiment  in  chronological  order  (e.g.,  in  BS,  depth 
a(675)  and  UML  was  14  and  8.8.  m  in  the  ‘shallow’  case  and  28  and  23.2  m  in  the  ‘deep’  case).  In 
general,  drastic  vertical  changes  of  a(675)  were  observed  in  the  vicinity  of  the  pycnocline  but  during 
the  FS  experiment  in  ES.  In  this  case,  it  is  suggested  that  phytoplankton  communities  were  not  actively 
growing  and  were  probably  sinking  as  part  of  a  post-bloom  stage  mainly  composed  df  senescent  cells 
(comm.  Pers.  Dr.  Percy  Donaghay).  Despite  these  differences,  spatial  patterns  of  ocean  color  above  the 
sea  surface  always  reflected  similar  modifications  when  optical  submarine  layering  became  deeper. 
Indeed,  the  skewness  of  R1  switched  from  positive  to  negative,  as  the  optical  submarine  layer  was 
placed  farther  from  the  sea  surface.  Also,  these  spatial  changes  in  ocean  color  can  be  attributed  to 
variability  of  lOPs  in  the  layer  above  the  main  ‘optieline’  due  to  the  relatively  shallow  penetration 
depth  of  passive  sensors  in  these  waters  (i.e.,  Zqd  always  smaller  than  depth  of  subsurface  a(675) 
layer). 

In  general,  skewness  of  a(675)  values  measured  above  the  subsurface  optical  layer  was  positively 
related  to  R1  skewness  but  in  Monterrey  Bay  experiments  (Fig.  3).  Likewise,  a(675)  above  and  within 
the  submarine  layer  tended  to  be  negatively  correlated  (Spearman  correlation  coefficient  up  to  -0.90  in 
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NS  of  ES)  but  in  MB  where  both  quantities  were  often  covarying  in  a  direct  way  (data  not  shown). 
Based  on  these  findings,  it  is  suggested  that  statistical  distribution  of  phytoplankton  cells  above  the 
main  ‘opticline’  was  linked  to  variability  of  lOPs  inside  the  submerged  layer  and  was  a  major  factor 
explaining  spatial  patterns  of  Rl.  The  discrepancy  found  in  Monterrey  Bay  may  be  attributed  to 
changes  on  phytoplankton  composition  (e.g.,  replacement  of  diatoms  by  dinoflagcllates,  Moline  et  al., 
2010).  However,  additional  experiments  arc  required  to  verify  this  hypothesis. 

For  each  study  area,  skewness  of  b  and  o  at  Rl  wavelengths  was  compared  between  NS  and  FS 
experiments  in  order  to  elucidate  which  optical  properties  and  wavelengths  were  responsible  of  Rl 
skewness  changes  due  to  depth  variation  of  a  single  submarine  optical  layer.  In  summary,  skewness  of 
Rl  appeared  to  be  mainly  connected  to  spatial  variability  of  total  scattering  and  presumably 
backscattcring  of  particulates  above  the  ‘opticline’.  These  optical  modifications  resembled  underwater 
light  field  perturbations  caused  by  internal  waves  (Weidemann  el  al.,  2001).  Lastly,  vj/Rl  changes  were 
mainly  driven  by  fe(440)/o(440)  rather  than  by  o(488)/6(488)  spatial  variability. 

Thermo-halinc  structure  as  inferred  from  vertical  distribution  of  seawater  density  was  usually  coupled 
to  depth  changes  of  biological  properties  such  as  a(675)  (Fig.  4).  Dekshenicks  et  al.  (2001)  found  that 
in  the  majority  of  cases  (i.e.,  >70%),  the  depth  distribution  of  optical  layers  is  mainly  determined  by 
the  water  column  stratification. 

Despite  this  overall  bio-physical  relationship,  decoupling  may  occur  (e.g.,  ES  in  Fig.  4)  when  optical 
constituents  do  not  behave  as  passive  tracers  (e.g.,  migrant  phytoplankton)  or  they  form  aggregates  that 
escape  the  pycnocline  barrier.  These  ecological  scenarios  may  introduce  a  large  uncertainties  when 
vertical  localization  of  submarine  optical  layers  is  based  on  hydrographic  profiles  (Zawada  et  al., 

2005).  Conversely,  vj/Rl  was  very  sensitive  to  vertical  changes  of  subsurface  optical  layers  for  a  broad 
range  of  vertical  mixing  conditions.  Thus,  this  study  supports  the  generalized  use  of  yRl  as  a  non- 
invasivc  optical  proxy  for  screening  relative  depth  (e.g.,  <20  m  versus  >20  m  depth)  of  thin  layers  in 
marine  waters  having  a  relatively  large  water  column  stability. 
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V  a(675) 


Figure  3.  Relationship  between  spatial  patterns  of  modeled  R1  and  field  measurements  of  a(675) 
above  the  main  ‘opticline’,  iff  is  the  skewness  of  the  empirical  or  theoretical  data  along  the  ROTV 
transect  \ff=0  corresponds  with  a  symmetrical  Gaussian  probability  distribution  (short-dash 
lines).  BS:  Black  Sea,  MB:  Monterrey  Bay,  ES:  East  Sound,  NS  and  FS  are  near-surface  and 

far-from-the  surface  case  studies,  respectively. 
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